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The sheer number of pulsars discovered by the SKA, in combination with the exceptional timing precision it 
can provide, will revolutionize the field of pulsar astrophysics. The SKA will provide a complete census of pulsars 
in both the Galaxy and in Galactic globular clusters that can be used to provide a detailed map of the electron 
density and magnetic fields, the dynamics of the systems, and their evolutionary histories. This complete census 
will provide examples of nearly every possible outcome of the evolution of massive stars, including the discovery 
of very exotic systems such as pulsar black-hole systems and sub-millisecond pulsars, if they exist. These exotic 
systems will allow unique tests of the strong field limit of relativistic gravity and the equation of state at extreme 
densities. Masses of pulsars and their binary companions — planets, white dwarfs, other neutron stars, and black 
holes — will be determined to ~ 1% for hundreds of objects. With the SKA we can discover and time highly-stable 
millisecond pulsars that comprise a pulsar-timing array for the detection of low-frequency gravitational waves. 
The SKA will also provide partial censuses of nearby galaxies through periodicity and single-pulse detections, 
yielding important information on the intergalactic medium. 



1. Introduction 

Neutron stars (NSs) are accessible to observa- 
tion as pulsars and thus provide our only means 
for probing the most extreme states of matter in 
the present-day Universe, which in turn will en- 
able a vast range of transforming science goals to 
be addressed, among which are 

• Strong-field tests of gravity including the 
Cosmic Censorship Conjecture and the no- 
hair theorem of BHs. 

• Detection of a cosmological gravitational 
wave background. 

• Mapping the complete structure of the 
Milky Way and revealing properties of the 
Galactic Center. 
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• Probing the intergalactic medium in new 
ways. 

• Identifying the equation-of-state of super- 
dense matter. 

• Quantifying the roles of magnetic fields and 
turbulence in core-collapse physics. 

• Understanding the superfluid interiors and 
relativistic magnetospheres of NS. 

• Unraveling the evolutionary and dynamical 
histories and properties of all Galactic glob- 
ular clusters. 

• Discoveries of extrasolar planets. 

Internal densities ~ ten times nuclear density 
have not existed since the Universe was about 1 
ms old. The actual state of matter in the core 
of a NS is presently not known. It may consist 
of de-confined quark matter or hyperonic matter 
produced in a phase transition that occurred dur- 
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ing or shortly after the core collapse of the pro- 
genitor star. Intermediate regions of the NS con- 
sist of neutrons and trace protons in, respectively, 
superfluid and super-conducting states achieved 
after the NS cooled to about 10 9 degrees. The 
outer regions include an ~ 1 km thick crust com- 
posed of iron-like nuclei surrounded by an ocean 
about 1 cm deep [66] ■ The magnetic field an- 
chored to the crust and extending to interstellar 
space is sufficiently strong that it elongates the 
atoms comprising the crust. The surface gravity, 
about 10 9 times that of the Earth's, is the largest 
of any object in the Universe subject to observa- 
tion, and corresponds to a gravitational redshift 
- 0.3. 

While NS gravity is strong, radio pulsars and 
probably most NSs are even more extreme elec- 
tromagnetically. By virtue of spin periods P ~ Is 
and magnetic fields B ~ 10 s — 10 14 Gauss, the 
electric force on a surface proton ~ 10 11 times 
larger than the gravitational force. Voltage drops 
~ 10 12 volts across the magnetosphere accelerate 
particles that can radiate across the entire electro- 
magnetic spectrum. This makes some pulsars vis- 
ible in every astronomical window |76j . Most of 
what is to be learned from pulsars requires obser- 
vations of their radio emission, often providing a 
unique source of information, and otherwise pro- 
viding information that complements multiwave- 
length studies, particularly at high energies. 

The feature that gives pulsars their name - 
pulsed radio emission — allows most NSs to be 
detected at levels well below the sky confusion 
limit and also provides the means for using pul- 
sars as physics laboratories. Coherent radio emis- 
sion is associated with the collimation of the flow 
of particles at the poles of the large-scale mag- 
netic field in combination with relativistic beam- 
ing. The spin-driven sweep of the beam across 
the line of sight then provides the distinctive pul- 
sation of the electromagnetic signal. 

The suitability for using pulsars as clocks de- 
pends on the regularity by which the spin evolves 
with time. Spin noise in some objects, which evi- 
dently reflects activity within the crust and su- 
perfluid, is large enough to mask many of the 
physical effects of interest that provide only sub- 
tle timing signatures. However, spin noise itself, 



including rapid spinups (glitches), is itself valu- 
able information on NS interiors. All pulsars are 
stable enough so that timing measurements can 
yield fundamental parameters such as the period 
vs. epoch, P(t) and its time derivative P, along 
with the dispersion measure DM. Objects with 
the narrowest pulses, the shortest periods, and 
the most stable rotation rates — millisecond pul- 
sars — yield the greatest opportunities for explor- 
ing relativistic gravity. Extrinsic gravitational ef- 
fects include perturbations of pulse arrival times 
from the direct influence of companion stars on 
space-time and from evolution of the orbit in the 
non-Newtonian gravity. The latter causes binary 
pulsars and their orbits to precess whilst their or- 
bit decays owing to loss of energy to gravitational 
radiation. Besides being sources of gravitational 
wave emission j^Z], pulsars also lend themselves 
as detectors of long-period gravitational waves 
that are cosmological in origin |19| . 

The population of isolated and binary pulsars 
is of great interest because their phase-space dis- 
tribution and overall numbers reflect the rate at 
which NSs are born in Type II supernova explo- 
sions, how the explosions themselves produce the 
runaway velocities of NSs, and how NSs are influ- 
enced by accretion in those rare binary systems 
that survive the explosions. Unlocking the vast 
population of active radio pulsars therefore allows 
us to study the star formation history and evolu- 
tion of massive stars, aspects of binary evolution, 
core collapse physics, as well as the movement 
of the high-velocity population of pulsars in the 
Galactic gravitational potential. 

With this enormous range of fundamental 
physics accessible through pulsar observations, 
the pulsar field has been extremely fruitful in the 
37 years since their discovery, as evidenced by 
the awarding of two Nobel prizes, one to the orig- 
inal discovery of pulsars , the other to the dis- 
covery of the first NS-NS binary pulsar that 
allowed inference of gravitational radiation in ac- 
cord with Einstein's General Theory of Relativity 
|75j . Nonetheless, the field has a great deal more 
to contribute to our knowledge of fundamental 
issues in physics and astrophysics. 

We describe some of the applications of pul- 
sars made possible by the Square Kilometer Array 
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project in the following Sections. In particular, 
we have identified Strong-Field Tests of Gravity 
Using Pulsars and Black Holes as a key science 
area for the SKA that is outlined in some detail 
in H3.ll and in particular in Kramer et al. (this 
volume) . 

To enable this programme of research outlined 
below, the SKA must have capabilities that al- 
low its enormous collecting area to be used in a 
variety of observing modes. Some of these imply 
usage as a huge, effective single dish, while others 
exploit the array aspects of the telescope. As a 
mantra for pulsar research, the SKA must allow 
us the following activities on pulsars: find them, 
time them and VLBI them. 

2. The Cosmic Census for Pulsars 

Pulsars are of great utility, no matter where we 
find them. So far, most known pulsars are Galac- 
tic, residing in or near the disk of the Galaxy or 
in globular clusters. A small number of pulsars 
is known in the Magellanic clouds; however radio 
pulsars have not yet been detected in more distant 
galaxies, although a few NS in accreting systems 
have been seen in M31 via their X-ray emission 
|38| . In the following we summarize why it is im- 
portant to discover more pulsars, both locally and 
in other galaxies, as well as in particular regions 
of the Milky Way. 

2.1. Galactic Census 

Why perform a full- Galactic census? The 
Galactic Census is essential in providing the lab- 
oratories for a wide range of pulsar applications 
that will be discussed in more detailed in the fol- 
lowing sections. The larger the number of pulsar 
detections, the more likely it is to find rare objects 
that provide the greatest opportunities as physics 
laboratories. These include binary pulsars with 
black hole companions (see H3.1J) : binary pulsars 
with orbital periods of hours or less that can be 
used for fundamental tests of relativistic gravity 
(see H3.1[) : MSPs that can be used as detectors of 
cosmological gravitational waves (see H3.3JI ; MSPs 
spinning faster than 1.5 ms, possibly as fast as 
0.5 ms, that probe the equation-of-state under 
extreme conditions (see H3.2I) ; hypervelocity pul- 
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Fi gure 1 . The P — P diagram for radio pulsars and 
magnetars. This is a scatter plot of pulse period and 
period derivative for pulsars in the current (Apr 2004) 
public ATNF catalogue. The circled points are pulsars in 
binary systems while objects in the upper right part of 
the diagram have derived surface magnetic fields exceed- 
ing Bq = 4.4 X 10 13 Gauss (see 813.41 . The "spin-up" line 
denotes the terminal period for objects moving to the left 
in the diagram owing to accretion that become "recycled" 
pulsars once accretion ceases. The "graveyard" marks the 
locations where radio pulsars emit less radiation or turn off 
entirely, possibly associated with the cessation of electron- 
positron (e^ ) pair-production cascades or with saturation 
of the radio luminosity at some fraction of the total spin- 
down energy loss rate. 
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sars with translational speeds in excess of 10 3 km 
s -1 , which probe both core-collapse physics and 
the gravitational potential of the Milky Way (see 
^13 . l(jpi : and objects with unusual spin properties, 
such as those showing discontinuities ( "glitches" ) 
and apparent precessional motions (both isolated 
pulsars showing 'free" precession and binary pul- 
sars showing geodetic precession) (see fl3.2|l . 

The second reason to perform a full Galactic 
census is that the large number of pulsars can be 
used to delineate the advanced stages of stellar 
evolution that lead to supernovae and compact 
objects. In particular, with as large a sample as 
possible, we can determine the branching ratios 
for the formation of canonical pulsars and mag- 
netars. We can also estimate the effective birth 
rates for MSPs and for those binary pulsars that 
are likely to coalesce on time scales short enough 
to be of interest as sources of periodic, chirped 
gravitational waves. The latter population can be 
directly compared to the results of gravitational 
wave detectors, which will have been operating 
for a sufficiently long time by the time the SKA 
is switched on. 

The third reason is that a maximal pulsar sam- 
ple can be used to probe and map the interstellar 
medium (ISM) in a nearly complete way. Measur- 
able propagation effects include dispersion, scat- 
tering, Faraday rotation, and HI absorption that 
provide, respectively, line-of-sight integrals of the 
free-electron density n e , of the fluctuating elec- 
tron density, 5n e , of the product B^n e , where B\\ 
is the LOS component of the interstellar magnetic 
field, and of the neutral hydrogen density. The 
resulting measures are DM, SM, RM and JVhi, 

DM = / dsn e , SM= / dsC 2 , (1) 
Jo Jo 

RM = / dsn c B lb N m = / dsn m , (2) 
Jo Jo 

involving quantities that can also be studied by 
other means, for instance by HI observations with 
the SKA. The determination of these observables 
for a large number of independent line-of-sights 
for pulsars will enable us to construct a complete 
map of the Galaxy (see UjTJJ- 

With a Milky Way birth rate that currently 



may be as large as 1CP 2 yr , about 10 8 NSs have 
been formed over the lifetime of the Galaxy, and 
probably many more because the star-formation 
was most likely higher in the past. Most NSs are 
inert, their radio emission having shut off long 
ago, and up to about half of them will have left 
the Galaxy owing to their large space velocities. 
Of the known pulsars (see Fig. , we can identify 
several subclasses: 

1. Canonical pulsars: These pulsars, like those 
first discovered, have present-day spin pe- 
riods ranging from tens of milliseconds to 
8 s and surface magnetic field strengths 
B ~ 10 12±1 G. They are often thought to 
be born with periods ~ 10 ms, though ev- 
idence suggests that some objects are born 
with periods longer than 0.1 s I n the 
standard picture of NS formation, all pul- 
sars start as canonical pulsars. In the P — P 
diagram of Fig. ^ most of these pulsars are 
located at P ~ Is and P ~ 10~ 15 . Young 
pulsars are especially important members of 
this class because they are associated with 
supernova remnants and often show copious 
numbers of glitches. 

2. Modestly recycled pulsars: are objects in bi- 
naries that survived a first SN explosion and 
subsequently accreted matter that spun-up 
the pulsar and reduced the effective dipolar 
component of the magnetic field. Accretion 
is terminated in these objects by a second 
supernova explosion that may or may not 
disrupt the binary. Those that survive are 
seen today as relativistic NS-NS binaries. 
Evolutionarily, it is possible that some sur- 
viving binaries include black-hole compan- 
ions. In the P — P diagram of Fig. these 
pulsars are typically located around P ~ 30 
ms and P ~ 10~ 18 . 

3. Millisecond pulsars (MSPs): objects in bi- 
naries that survive the first SN explosion 
and in which the companion object eventu- 
ally evolves into a white dwarf. The long, 
preceding accretion phase spins the pulsar 
up to millisecond periods while attenuating 
the (apparent) dipolar field component to 
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10 8 — 10 9 G. The consequent small spin- 
down rates seem to underly the high timing 
precision of these objects and imply spin- 
down time scales that exceed a Hubble time 
in some cases. In the P — P diagram of 
Fig. ^ these pulsars are typically located 
around P ~ 5 ms and P ~ 10~ 20 . Evolu- 
tionary scenarios that produce recycled pul- 
sars and MSPs are discussed in UJ. 

4. Strong-magnetic-field pulsars: Recently dis- 
covered radio pulsars have inferred fields 
>10 14 G 852, rivalling those inferred for 
"magnetar" objects identified through their 
X-and-7 radiation that seems to derive from 
non-rotational sources of energy. The rela- 
tionship between magnetars and these high- 
field radio emitting pulsars, whose radiation 
derives solely from spin energy, is not yet 
known. In the P—P diagram of Fig.^thcsc 
radio pulsars are typically located around 
P ~ 5s and P ~ KT 13 . 

We envision a full-Galactic census of radio pul- 
sars that aims to detect at least half of the ac- 
tive radio pulsars that are beamed toward Earth 
(Figure |2J. The typical lifetime of canonical pul- 
sars ~ 10 Myr, where we define lifetime as the 
duration of the radio-emitting phase that, for ob- 
jects with B ~ 10 12 G, is the time needed for a 
rapidly rotating object to reach the "death band" 
on the right-hand side of the P — P diagram. At 
long periods short-ward of the death band, where 
pulsars spend most of their detectable lifetimes, 
the beaming fraction ~ 20% (e.g. [22), so the 
fiducial birth rate implies ~ 2 x 10 4 detectable 
pulsars in the Galaxy. Non-canonical classes of 
pulsars add to these numbers only negligibly be- 
cause their effective birth rates are smaller by a 
factor 10- 4 -10~ 3 . 

2.2. Galactic Center 

The Galactic Center (GC) is an especially tan- 
talising but exceedingly difficult region to search 
for pulsars. In many respects, the GC is similar 
to an especially large globular cluster with regard 
to the density of stars (cf . H2.3(l . It differs in that 
molecular material and, hence, star formation are 
both much more prominent in the GC than in 
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Figure 2. Simulated Galactic pulsar population discov- 
ered in a SKA survey of the enitre sky. The ~ 20, 000 pul- 
sars shown together with the spiral arms structure. The 
Galactic Center is located at the origin while the Sun is 
at (0.0,8.5) kpc. 



globulars. Additionally, the ~3x 1O 6 M black 
hole 64 that underlies the compact radio source, 
Sgr A*, perturbs space time significantly and is 
fed episodically by inspiraling gas and stars. 

Why find pulsars in the GC? Any pulsars de- 
tected in or beyond the GC that are viewed 
through the region are potentially unique probes 
of the gas, magnetic field, and space-time of the 
GC region. First of all, pulsars in the GC can 
be used to probe the magnetoionic medium along 
the line-of-sight, including possible detection of 
the inner scale for electron density turbulence (see 
i j3.8|) . Second, the initial mass function and over- 
all stellar evolution in the GC is likely to be very 
different from the rest of the Galaxy. The num- 
ber and ages of pulsars and their binary member- 
ship will provide clues about these areas (see 
and ij3. 1U|) . Finally, the large stellar density offers 
the possibility of finding pulsars with stellar black 
hole companions, allowing unprecedented tests of 
gravitational theories in the strong field limit and 
the study of black hole properties. Similar stud- 
ies will be possible for Sgr A*, the supermassive 
black hole in the center, if pulsars in compact or- 
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bits around the black hole are found (see H3.1fl . 
Pulse timing measurements may provide the pos- 
sibility of measuring the spin of the GC's central 
black hole. 

At present, none of the known pulsars is within 
or beyond the GC. The primary reason is that an 
especially intense scattering region lies between 
us and the GC at close proximity (~ 100 pc) to 
the GC ■ The large scattering has been known 
since shortly after Sgr A* itself was discovered 30 
years ago and it has been probed through angular 
diameter measurements of OH/IR masers in the 
GC region and through surveys for background 
AGNs. The implication for pulsars is that, at 
a standard search frequency of 1.4 GHz, a pulse 
emanating from the location of Sgr A* would be 
broadened to ~ 300 s owing to multi-path prop- 
agation. 2 

The pulse broadening time Ta scales with the 
observed angular diameter 9d as = W T D6i/2c, 
where D is the source distance and W T is a geo- 
metrical weighting factor that takes into account 
the location of the scattering region relative to 
the source and observer. For the GC, the weight- 
ing factor is very large because the region is much 
nearer the GC than the observer. For scattering 
by fluctuations in the electron density, 6d oc v~ 2 
and thus the broadening time is a very strong 
function of frequency, oc v~ A . 

To combat pulse broadening, observations at 
higher frequencies are therefore needed that ex- 
ploit the strong v~ A dependence. For example, 
10 GHz observations yield ~ 0.11s broadening, 
small enough to allow detection of longer period 
pulsars. However, pulsar spectra are power-law in 
form, often with steep dependences tx v~ a with 
a ranging from to 3 and (a) w 1.5. While ongo- 
ing searches may yield detections of a few pulsars 
in the GC, some of which may provide impor- 
tant probes of the region, the sensitivity of the 
SKA is clearly needed at high frequencies to de- 
tect a meaningful sample of pulsars in the GC, as 
demonstrated in Figure The figure shows the 
radio luminosity (defined as L p = flux density x 
distance 2 ) at 1.4 GHz plotted against spin period 




2 Pulsars located closer to us than Sgr A* but still behind 
the screen are scattered less than this while pulsars beyond 
the GC viewed through the screen are more scattered. 
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Figure 3. Detectability of pulsars in the Galactic Cen- 
ter and in other galaxies is shown in this plot of "pseudo 
luminosity" L p = ShqoD 2 vs. spin period for pulsars in 
the current (Apr 2004) public ATNF catalogue; S1400 is 
the period-averaged flux density at 1.4 GHz. Given duty 
cycles for pulsars ranging from 0.01 to 0.7, the peak flux 
density can be a large multiple of the period-averaged flux 
density. Solid and dashed lines mark the minimum de- 
tectable spin period P m i n for pulsars at the Galactic Cen- 
ter that are detectable with the GBT at 9 GHz and the 
SKA at 9 and 15 GHz. Scattering is assumed to be in a 
filled region centered on Sgr A* with 1/e cylindrical radius 
of 0.15 kpc, as included in the NE2001 model. For pulsars 
at the location of Sgr A*, the dispersion and scattering 
measures are DMqc = 1577 pc cm -3 and SMqc = 10 7 
kpc m" 20 / 3 . The curves are based on harmonic sums 
of the Fourier power spectra of dedispersed time series. 
They correspond to surveys at these frequencies but have 
all been scaled to the equivalent sensitivity at 1.4 GHz 
assuming that the pulsar flux densities scale as oc u~ 2 . A 
bandwidth of 800 MHz was assumed in all cases. Inte- 
gration times of 1.5 hr and 10 4 s were used for the GBT 
and the SKA, respectively. In addition, the SKA curves 
assume that the full sensitivity of the SKA is available. 
Dotted horizontal lines show the minimum values for L p 
needed for detection at 1.5 GHz in periodicity searches of 
pulsars at distances of 1 Mpc and 5 Mpc. The dotted lines 
assume detectability is independent of spin period, which 
will not be the case if there is significant scattering along 
the line of sight. 
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along with detection threshold lines for the Green 
Bank Telescope (GBT) and for the SKA. 
The maximum detectable distance is 

D = f Lp V 2 N 1/4 

"max — I I Jv h ) 

\ ^mini / 

where S mini — mS sys /V2BT is the single- 
harmonic threshold with m = the threshold 
signal-to-noise ratio, and is the number of har- 
monics detected in the search power spectrum. 
The minimum luminosity for m = 10 is thus 

tpmlI = 8 T.5„w(^) 2 (f)"' /2 . 

For short period pulsars, the number of harmon- 
ics detected will be less than the 16 number as- 
sumed in this equation. 

Also shown in Figure OH are detection curves 
for the minimum detectable period P m in vs. 
L p (1.4GHz). These are shown for pulsars at 
the Galactic Center that are detectable with the 
Green Bank Telescope (GBT) at 9 GHz and with 
the SKA at 9 GHz and 15 GHz. These curves 
were calculated for pulsars near the location of 
Sgr A* and take into account pulse broadening 
from scattering and the spectral dependence of 
the pulsar flux density using a conservative value 
for spectral index, a = 2. The detectability was 
thus calculated for surveys at these frequencies 
but has been scaled to the equivalent sensitivity 
at 1.4 GHz. The calculations assume that the 
full sensitivity of the SKA is available. If only a 
fraction f c can be used in a blind survey using a 
"core" array, then the curves will move upward 
by an amount log/ c . 

2.3. Globular Clusters 

Globular Clusters hold vast reservoirs of MSPs, 
which are formed at a rate per unit mass which is 
at least an order of magnitude higher than in the 
Galactic disk. The reason for this overabundance, 
which may be even higher than presently seen due 
to the escape of high velocity objects, is thought 
to be the formation of binaries via two- and three- 
body encounters in the high density environments 
of Globular Clusters. 

Why perform a full census of globular clusters? 
The greatly increased likelihood that many of 



these MSPs in globular clusters will have under- 
gone some form of dynamical interaction means 
that the chances of finding exotic binaries, such 
as the long sought-after MSP-black hole system, 
are perhaps highest in globular clusters (see H3.1|) . 
Furthermore the pulsars in each globular cluster 
will provide us with exceptional probes of the his- 
tory and evolution as well as its present prop- 
erties, including the dynamics, gas content, ac- 
curate distances and proper motion (see ^3.91) . 
Moreover, an open question is whether globular 
clusters contain massive black holes in their cen- 
ters |54I31| . or possibly even binary black holes 
|12| . Pulsar timing with the precision achievable 
with the SKA will provide us with a tool to reveal 
whether such systems are present. This can be 
achieved either by probing the inner most regions 
of the cluster to reveal velocities or accelerations 
expected due to the presence of a black hole, or by 
using the pulsars in the cluster as an in situ grav- 
itational wave detector sensitive to the presence 
of binary black holes (see ^3.3[) . 

At present there are 76 pulsars known in 23 
globular clusters. Not all clusters are equally rich, 
and it is not clear what determines the number 
of active pulsars in a given cluster. Based on our 
knowledge of the luminosity function of the MSPs 
in the Galactic disk combined with the continuum 
radio emission and the numbers of MSP-like X- 
ray sources associated with globular clusters, the 
total population is likely to be at least two orders 
of magnitude more than this. 

The remaining pulsars in these globular clus- 
ters are difficult to find with present instruments 
due to a combination of the intrinsic low luminos- 
ity of MSPs and the typically large distances to 
globular clusters. Furthermore, most pulsars are 
in compact binary orbits. Given the presently 
long integration times needed to detect the pul- 
sars, the pulse signal is often smeared out due to 
the Doppler effect. Recent developments of so- 
phisticated techniques to correct for this smear- 
ing have resulted in an increase in the number of 
systems known, but blind searches are often too 
computationally expensive and the vast majority 
of systems are still too weak for detection. 

Targeted SKA surveys of all Galactic globular 
clusters will enable us to detect all of the pulsars 
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contained therein which are beamed in our di- 
rection. This complete census will uncover large 
numbers of pulsars with which to study their for- 
mation, evolution, spin parameters, binary nature 
and emission properties. 

2.4. External Galaxies 

Pulsars beyond the disk of the Milky Way are 
known only in globular clusters and in the Mag- 
ellanic Clouds owing to their intrinsic faintness. 
With the SKA, galaxies in the local group, includ- 
ing M31 and M33, are within reach using peri- 
odicity searches while giant pulses like those seen 
from the Crab pulsar can be detected from galax- 
ies out to at least 5 Mpc. 

What is the importance of detecting pulsars in 
other Galaxies? Pulsars likely to be detected 
will be young pulsars with high luminosities that 
can be correlated with catalogs of supernova 
remnants and will yield estimates of the star- 
formation rate and the branching ratio for su- 
pernovae to form spin-driven pulsars as opposed 
to magnetars and blackholes. Extragalactic pul- 
sars will also provide information about the mag- 
netoionic media along the line of sight through 
determination of DM, SM, and RM. Unambigu- 
ous study of the intergalactic medium in the local 
group requires removal of contributions to these 
measures from the foreground gas in the Galaxy 
and gas in the host galaxy. The more pulsars de- 
tected in a galaxy, the more robust this removal 
will be. 

Extragalactic pulsars can be found through 
blind surveys for both periodic sources and in- 
dividual giant pulses. Additional successes will 
follow from targeted surveys of individual super- 
nova remnants in the nearest galaxies. In M.ll we 
discuss the requirements on sensitivity and field 
of view (FOV). 

Figure includes detection lines for surveys at 
1.4 GHz for pulsars at 1 Mpc and 5 Mpc assuming 
that the full SKA sensitivity is available. Full 
sensitivity would apply to targeted searches of, for 
instance, supernova remnants in nearby galaxies. 
However, ~full-FOV blind surveys will provide 
only a fraction f c corresponding to a core array. 
With full sensitivity, a reasonable fraction of the 
luminosity function can be sampled out to 1 Mpc 



and a few pulsars can be detected to 5 Mpc. 

Giant pulses from the Crab pulsar serve as 
a useful prototype for estimating detection of 
strong pulses from nearby galaxies. The strongest 
pulse observed at 0.43 GHz in one hour has 
5/A max = 10 4 - even with the system noise 
dominated by the Crab Nebula. For objects in 
other galaxies, the system noise is dominated by 
non-nebular contributions, implying that the S/N 
in this case would have increased by a factor of 
about 300. We can estimate the maximum dis- 
tance of detection at a specified signal-to-noise 
ratio, (S/N) det as 



1.6 Mpc 



y/(S/N) dot /5 \A Arecih J 



1/2 



(3) 



where / c ^4ska is the SKA's collecting area that 
can be used for a giant pulse survey and A&recibo 
is Arecibo's effective area at 0.43 GHz. The us- 
able area for the SKA in a blind survey will 
be limited by what fraction f c of the anten- 
nas are directly connected to a central correla- 
tor/beamformer. For f c = 1, AsKA/^Arecibo ~ 
20, the standard one-hour pulse seen at Arecibo 
could be detected out to 7.3 Mpc. 

3. Fundamental Physics &; Astrophysics 

Having discovered a large sample of pulsars in 
the census of the Milky Way, the Galactic Cen- 
tre, Globular Cluster and external galaxies, a vast 
range of fundamental problems in physics and as- 
trophysics can be studied. We highlight some of 
those in the following. 

3.1. Tests of Theories of Gravity 

The theory of General Relativity (GR) has to 
date passed all observational tests with flying 
colours. Nevertheless, one of the most funda- 
mental questions remaining is whether Einstein's 
theory is the last word in our understanding of 
gravity or not. Solar system tests of GR are 
made under weak-field conditions, illuminated by 
the small orbital velocities and gravitational self- 
energies involved (e.g. the self-energies for the 
Sun and the Earth, expressed in units of their 
rest-mass energy, are es U n ~ — 10 -6 and CEarth ~ 
— 10~ 10 , whilst for a pulsar and black hole (BH) 
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we find epsr ~ —0.2 and €bh = —0.5). There- 
fore, none of these tests of gravitational theories 
can be considered to be complete without prob- 
ing the strong-field limit, which is done best and 
with amazing precision using binary radio pul- 
sars |18|. However, even the existing binary pul- 
sar tests only begin to approach the strong-field 
regime and the discovery of more extreme binary 
systems is required. Indeed, the SKA is the only 
instrument which promises to probe the strong- 
field limit via the discovery and timing of such 
systems, in particular that of a pulsar orbiting a 
BH. The promises provided by the discovery of 
such particular systems are discussed in detail in 
a companion key science article (Kramer et al., 
this volume) , where we detail how pulsars can be 
used to test the "Cosmic Censorship Conjecture" 
and the "No-hair" theorem for the description of 
BHs. Here, we concentrate on general aspects of 
tests of gravitational theories using pulsars timed 
with the SKA. 

Since NSs are very compact massive objects, 
double neutron star (DNS) binaries can be con- 
sidered as ideal point sources. Finding and timing 
DNSs in tight binary orbits - ideally close to co- 
alescence and thus emitting strong gravitational 
waves - provide stringent tests of theories of grav- 
ity in the strong-gravitational-field limit. Tests 
can be performed when a number of relativistic 
corrections to the Keplerian description of the or- 
bit, the so-called "post-Keplerian" (PK) param- 
eters, can be measured. For point masses with 
negligible spin contributions, the PK parameters 
in each theory should only be functions of the a 
priori unknown NS masses and the well measur- 
able Keplerian parameters. With the two masses 
as the only free parameters, the measurement of 
three or more PK parameters over-constrains the 
system, and thereby provides a test ground for 
theories of gravity [UJ . In a theory that describes 
the binary system correctly, the PK parameters 
produce theory-dependent lines in a mass-mass 
diagram, which compares the masses of the two 
neutron stars, that all intersect in a single point. 

The best example for such tests is cur- 
rently given by the double pulsar system PSR 
J0737— 3039 where five PK parameters are avail- 
able for tests, in addition to the theory indepen- 



dent mass ratio of the two NSs For the rel- 
ativistic binary systems that will be discovered 
and timed with the SKA, one would be able to 
routinely measure five or more PK parameters, 
severely over-constraining the system. In partic- 
ular, more double pulsar systems will be found, 
and PK parameters that are currently impossible 
(or extremely difficult) to measure, such as those 
arising from aberration effects and geodetic pre- 
cession, would become accessible with the SKA, 
providing also full 3-D information about the ori- 
entation of these binary systems. 

Most importantly, current PK parameters are 
only measured to the lowest post-Newtonian ap- 
proximation. The timing precision achievable 
with the SKA means that higher-order correc- 
tions are likely to become important, demand- 
ing the development of a timing formula that is 
accurate to at least second post-Newtonian or- 
der. For instance, corrections would need to ac- 
count for the assumptions currently made in the 
computation of the Shapiro delay (see H3.2|) that 
gravitational potentials are static and weak every- 
where |83l43l44j . In addition, other effects such 
as those related to light bending and its conse- 
quences will become important where an addi- 
tional signal would be superposed on the Shapiro 
delay as a typically much weaker signal, which 
arises due to a modulation of the pulsars' rota- 
tional phase by the effect of gravitational deflec- 
tion of the light in the field of the pulsar's com- 
panion 20 . 

Like the DNS binaries, most NS-white dwarf 
(NS-WD) binaries can also be considered as pairs 
of point masses. Significant measurements of PK 
parameters can only be obtained in a few cases 
so far, since relativistic effects for NS-WDs are 
generally much smaller than for DNS binaries 
(e.g. [EHED- However, with the SKA, PK pa- 
rameters can be determined for many more sys- 
tems. Such measurements can be similarly valu- 
able since NS-WD systems test different aspects 
of gravitational theories. For instance, as shown 
recently [23], the tests provided by the PSR- WD 
system J1141— 6545 are more constraining for the 
class of tensor-scalar theories than those tests 
provided by the double pulsar. The reason is 
that unlike GR, some alternative theories of grav- 
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ity, such as tensor-scalar theories, predict effects 
that depend strongly on the difference between 
the gravitational self-energy of two orbiting bod- 
ies, which is large in NS-WD systems. Using such 
systems, one is able to probe possible violation 
of conservation of momentum, equivalence princi- 
ples and expected Lorentz- and positional invari- 
ances jTSj- A manifestation of such effects would, 
for instance, be measurable in the detection of 
gravitational dipole radiation where as GR only 
predicts quadrupole emission. 

The computer power available when the SKA 
comes online will enable us to do much more so- 
phisticated searches in parameter space, includ- 
ing full searches in acceleration due to binary mo- 
tion, than possible today, and the SKA sensitiv- 
ity allows much shorter integration times, so that 
searches for compact binary pulsars will no longer 
be limited. Hence, the combination of SKA sen- 
sitivity and computing power means that the dis- 
covery rate for relativistic binaries is certain to 
increase beyond the number of at least a hundred 
compact binary systems that we can expect from 
an extrapolation of the present numbers. 

In summary, tests of gravity affordable with 
the SKA will not simply be a continuation of the 
present tests at higher precision levels, but the 
better sensitivity and timing precision will en- 
able us to perform genuinely new tests. These 
prospects arise from the assumption that the tim- 
ing precision can be improved by two or even 
three orders of magnitude over the current stan- 
dards. However, as we discuss in some detail in 
this requires not only an increase in raw sen- 
sitivity but also the ability to correct for propa- 
gation effects from intervening plasmas and suf- 
ficient polarization purity and (quasi-) simulta- 
neous multi-frequency observations. Moreover, 
the intrinsic phase-jitter of pulsars will prevent 
the timing precision of some MSPs from reaching 
the theoretical limit given by radiometer noise, so 
that the application of correction schemes needs 
to be considered on a case-by-case basis. How- 
ever, the current timing precision of a number of 
MSPs appears to be limited by telescope sensi- 
tivity, as opposed to systematic effects (e.g. 
suggesting that much improved timing precision 
can be achieved for a few objects and that new, 



unique tests of relativistic gravity will be enabled 
through greater telescope sensitivity. 

It is important to note that many of these ex- 
periments require corrections of the measured pa- 
rameter values for kinematic effects. For instance, 
the precision of GR tests achievable with PSR 
B1913+16 is limited by the accuracy to which 
the pulsar's motion and acceleration in the grav- 
itational potential of the Galaxy is known |82|. 
Observed values for parameters like the orbital 
decay rate, Pb/Pb, or changes in the semi- major 
axis, x/x, are affected by a kinematic Doppler 
term given by |17| 

- -p: = - K ■ (a PS R - a SS Bj H -j , (4) 

D c cd 

where Kq is a vector from the Solar System 
Barycentre (SSB) towards the binary pulsar, 
apsR aud asSB are the Galactic accelerations at 
the location of the binary system and the SSB. 
The last term including the transverse velocity, 
Vr, and distance, d, to the pulsar is known as sec- 
ular acceleration or "Shklovskii term" [70]. Cor- 
recting for this term to derive the intrinsic val- 
ues obviously requires precise astrometric infor- 
mation like proper motion and distances which 
can be derived using the VLBI capabilities of the 
SKA. In order to get precise distances for pulsars 
across the Galaxy a final positional accuracy of 
0.1 mas is needed at 5 GHz. 

3.2. Structure of Neutron Stars &; 
Equation-of-State 

The internal structure of NSs is complex. This 
question was already addressed by Oppenheimer 
& Volkov (1939) long before the discovery of pul- 
sars. The structure depends sensitively on the 
equation of state, i.e. the relationship between 
density and pressure. Knowledge of the exact 
equation of state would allow us to deduce most 
of the NS's physical properties, most notably its 
radius for a given mass. With the SKA we can 
address this question from various angles. 

Firstly, the study of the equation of state will 
benefit from the vastly improved statistics of ob- 
served pulse periods. Periods of known pulsars 
range from 1.5 ms to 8.5 s. The discovery of 
even smaller rotational periods will provide sig- 
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nificant insight into the properties and stiffness 
of the ultra-dense liquid interior of NSs. Even 
the case of a "null result" , that no period shorter 
than 1.5 ms will be found, will require a theoret- 
ical explanation, as it requires the existence of a 
limiting period not far away from the currently 
observed value but larger than our current the- 
oretical limits. At the other extreme end, the 
discovery of very long period pulsars, isolated or 
in binary systems, will establish the connection 
of radio pulsars to magnetars, Soft-Gamma ray 
Repeaters (SGRs) and Anomalous X-ray Pulsars 
(AXPs), clarifying as to whether these are dis- 
tinct classes of NSs or simply different evolution- 
ary stages or end products (see also H3.4fl . 

Secondly, observations of relativistic effects in 
binary pulsars allows precise determinations of 
the masses of the pulsar and its companion. This 
is possible when either two or more PK param- 
eters can be determined (see H3.ll) . or when a 
Shapiro delay can be measured. The Shapiro de- 
lay in the arrival time of a pulse arises from the 
curvature of space-time caused by the presence of 
a companion star. Measuring the Shapiro time 
delay is therefore the most straightforward way 
to obtain a direct measurement of the compan- 
ion mass and, hence, of the mass of the active 
pulsar from knowledge of the total system mass 
(e.g. by measuring a relativistic precession of the 
orbit). Currently, the masses of about 20 neutron 
stars can be determined (Fig.^J. With hundreds 
of NS mass measurements available, the range of 
observed masses can be studied and related to 
the spin periods, promising to reliably establish 
the maximum possible mass of observable NSs. 
Moreover, the observed values can be related to 
the evolutionary history of the binary systems, 
studying the amount of matter accreted during a 
spin-up process as an X-ray binary. Given that 
some millisecond pulsars are very old, one may 
also be able to put limits on the evolution of NS 
masses with age of the universe, so that variations 
in the gravitational constant can be detected or 
constrained [ZZ|- It is also worth pointing out 
that, in addition to precisely measuring the dis- 
tribution of masses of NSs, the same observational 
techniques can be used to obtain mass measure- 
ments of other compact pulsar companions. The 



current timing precision usually prevents this, but 
with the SKA it would be possible to perform 
such measurements also for WD companions. fTQ 



Neutron star mass (Solar Masses) 

Figure 4. Mass measurements of neutron stars obtained 
from the timing of binary pulsars with corresponding un- 
certainties shown. The canonical value of 1.4Mq is in- 
dicated by the vertical line. While mildly recycled and 
young pulsars are broadly consistent with this value, mil- 
lisecond pulsars born in the process of a prolonged mass 
transfer appear to have somewhat larger masses. Figure 
compiled by Stairs (2004). 



Thirdly, the moment of inertia of NSs can be 
determined when high precision timing observa- 
tions allow the detection of relativistic spin-orbit 
coupling in DNSs. This effect would be visible 
most clearly as a contribution to the observed 
periastron advance in a secular 0j and periodic 
fashion [55] . The size of the contribution depends 
on the pulsars' moments of inertia, which there- 
fore can be determined if two other PK parame- 
ters can be determined with sufficient accuracy so 
that the effect of spin-orbit coupling can be iso- 
lated [23 ■ The SKA timing precision will be cru- 
cial in pushing the measurement accuracy to the 
limits required to see this effect of second Post- 
Newtonian (i.e. (v/c) A ) order (cf. jEH) - 



Double Neutron Star 



Recycled Pulsar-WiD Systems 



J0737-3039A (a) 
J0737 3039B (a) 
J1518-4904 (b) 
J1518+4904eom (b) 
B1534+12 (e) 
B1534+12com (c) 
J1B11-1736 (d) 
J1811-1736com (d) 
B1913+16 (e) 
B1913 ■ 16com (e) 
B2127+1 1C (f) 
B2127+llCcom (f) 



J0045-7319 (g) 
J1141-6545 (h) 
B2303 + 46 (i) 



J0437-4715 (j) 
J0621+1002 (k) 
J0751 ■ 1807 (1) 
J1012+5307 (m) 
J1713+0747 (n) 
B1802-07 (o) 
B1855+09 (p) 
J2019+2425 (q) 
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Finally, observations of rotational instabilities 
in pulsars, which include rapid spin-ups followed 
by slow relaxation — glitches — and stochastic 
variations generically called timing noise, probe 
the interior of NSs in a unique way that can be 
considered "neutron star seismology" . Observa- 
tions of glitches and their recovery provide strong 
evidence for the existence of a fluid component in- 
side the solid outer crust of the NS, providing the 
basic picture of NS structure provided in the in- 
troduction. This picture is based on only a small 
number of glitches observed in even fewer pulsars 
(e.g. [SHI)- The reason is that glitches - mostly 
observed in young pulsars - are rare, so that their 
timely detection is important to follow the reveal- 
ing relaxation process. With the SKA a large 
number of pulsars would be monitored simulta- 
neously, so that glitches can be observed as they 
happen. The multi-beam capabilities of the SKA 
are essential for this study. 

The results obtained from the analysis of the 
detected glitches can be compared to the obser- 
vations of pulsars that are found to be freely 
precessing, since the time-scales for precession 
are probably related to the coupling-strength be- 
tween the liquid interior and the solid crust. Cur- 
rently, convincing evidence for free precession is 
found for only a few pulsars but SKA tim- 
ing observations of all pulsars discovered in the 
cosmic census will produce further examples. 

Timing noise constitutes stochastic variations 
in pulse phase residuals (e.g. after removing spin- 
down polynomial components and orbital effects) 
that show non-stationary statistics and is thus of- 
ten described as having a "red" power spectrum. 
It is radio-frequency independent (after any DM 
and scattering variations from the ISM are sub- 
tracted) and appears to be caused by torque vari- 
ations acting on the NS crust. It is possible that 
some or all timing noise is associated with free 
precession induced by crust quakes that cause 
misalignment of the angular momentum vector 
and the spin axis. This conjecture can be tested 
only by identifying more objects that show pulse 
shape variations caused by wobble of the radio 
beam that accompanies pulse phase variations. 



3.3. Pulsars as sources and detectors of 
gravitational waves 

Pulsars are not only detectable across the 
whole electromagnetic spectrum, but rotating 
and coalescing NSs are also among the sources 
expected to be detected first with gravitational 
wave (GW) detectors. 

The recent discovery of the J0737— 3039 sys- 
tem suggests that Advanced LIGO, operating on 
timescales similar to those for the SKA, will de- 
tect a few coalescing binary pulsars each day [7| 
I39| . Moreover, it is anticipated that J0737— 3039 
can act as a calibrating signal for the space in- 
terferometer LISA. The SKA will discover many 
more such compact relativistic binaries and will 
hence provide key information for a location and 
study of these objects in a non-photonic window. 

Moreover, deviations from a purely spherical 
shape will cause spinning NSs to emit continu- 
ous GW emission that may even be detected with 
current generations of GW telescopes like LIGO, 
GEO600 or VIRGO. The faintness of the signal 
requires knowledge of the precise positions and 
spin-frequencies of the pulsars. With the major- 
ity of Galactic pulsars potentially being discov- 
ered with the SKA (see M2.1|l . sensitive targeted 
GW searches are possible for a very large number 
of pulsars, so that the actual shape of a NS can 
be probed. 

As a result of the cosmic census, the SKA will 
also produce a dense array of millisecond pulsars 
across the sky. Being timed to very high pre- 
cision (<100 ns), these act as multiple arms of a 
cosmic GW detector when a passing gravitational 
wave perturbs the metric and hence affects the 
pulse travel time and the measured arrival time 
at Earth |l9l2SI60| . This "device" with the SKA 
at its heart - the so-called "Pulsar Timing Ar- 
ray" (PTA) - will be sensitive to GWs to frequen- 
cies of 1/observing time, hence ^nHz. The PTA 
thereby complements the much higher frequen- 
cies accessible to Advanced LIGO (~100Hz) and 
LISA (~mHz), and the extremely low frequen- 
cies probed by polarization studies of the Cosmic 
Microwave Background (~ 10~ 18 Hz). 

The measurement precision and accuracy of the 
pulsar clock is not sufficient to detect the gravi- 
tational radiation of stellar-mass binaries by the 
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means of a PTA. However, super-massive black 
hole binaries in nearby galaxies with orbital pe- 
riods of a few years would produce a periodicity 
in pulsar arrival times with an amplitude of the 
order of 10 ns to 1 /is and in a frequency range 
that is detectable with SKA timing |60I49| . The 
light-travel time delay between the Earth and the 
timed pulsars can, depending on the geometry, 
enable us to observe a given massive binary sys- 
tem at two different epochs simultaneously. A 
slow decay of the binary orbit hence results in 
both low and high frequency components of the 
timing residuals, whereas the difference in the 
frequencies of these components will depend on 
the orbital decay rate. As pulsar timing is more 
sensitive to lower frequencies, the highest ampli- 
tude oscillations in the timing residuals will be 
due to the delayed component. This effect cor- 
responds to the three-pulse response occurring in 
spacecraft Doppler tracking experiments and the 
multi-pulse response from time-delay interferom- 
ctry used with LISA 36 . 

The short lifetime of such massive binaries re- 
duces the chances of detecting such systems. In 
contrast, the SKA can detect the signal of a 
stochastic background of gravitational wave emis- 
sion produced by a large number of unresolved 
independent and uncorrelated events. Another 
contribution to the GW background, but with a 
different spectrum, is expected in some cosmo- 
logical string theories. The SKA is crucial in an- 
swering the question about the existence, nature 
and composition of a stochastic GW background. 
The construction of a PTA to detect such signals 
is one of the key-science projects of the SKA and 
is presented in more detail in the contribution by 
Kramer, this volume. 

3.4. Relativistic Plasma Physics 

The plasma physics occurring under extreme 
conditions of very high densities, super-strong 
magnetic fields and in the rather complex envi- 
ronment known as the pulsar magneto sphere are 
only poorly understood. This situation remains 
after almost 40 years of pulsar research and is in 
spite of the large number of studies and observa- 
tions devoted to the identification of the relevant 
physical processes (e.g. 123)- Solving this "pul- 



sar problem" however offers deep insight into rela- 
tivistic plasma physics, and the SKA may provide 
us finally with the observational facts to under- 
stand pair production processes, the creation of 
radio and high-energy emission and the overall 
structure of the pulsar magnetosphere. 

Radio Fluctuations: The size of the region 
of electromagnetic activity around active pul- 
sars, the magnetosphere, is essentially the light- 
cylinder radius, 

^ = £ = ^4.77x10^^). (5) 

Radio emission altitudes are much smaller than 
i?LC and radio emission is expected to originate 
as close to the NS surface as 10-100 stellar radii, 
produced by a relativistic pair plasma that is cre- 
ated in the polar gap region |46|. 

While pulsars are renowned for their stable av- 
erage pulse profiles formed by the summing to- 
gether of many pulses, single pulses provide us 
with the rich detail that reflects the physics of 
the generation and propagation of radio emis- 
sion. The initially chaotic appearance of single 
pulses can often be resolved to reveal a rich phe- 
nomenology that includes: drifting sub-pulses, 
nulling, microstructure and giant pulses which 
correspond to the building blocks of the observed 
pulsar radiation (e.g. |79I32| ). Key tests of all 
emission mechanisms and the role of the mag- 
netosphere therefore come from determining the 
instantaneous time duration and evolution, the 
frequency bandwidth and the polarisation prop- 
erties of these phenomena. However, our under- 
standing of how this phenomenology reflects the 
underlying physics has been restricted because ra- 
dio emission from pulsars is weak and presently 
there are only a limited number of sources in 
which single pulses can be studied (e.g. The 
unprecedented leap forward in sensitivity com- 
bined with sky coverage provided by the SKA will 
revolutionise the study of single pulses. To com- 
plete this study, instantaneous coverage of a very 
wide band (e.g. 0.1 to 5 GHz or higher) is de- 
sireable so that individual pulses can be tracked. 
This specification surpasses the strawman 20% 
bandwidth specification for the SKA but would 
allow new approaches to the understanding of 
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pulsar magnetospheres. 

Relationship to Magnetars and Other Objects: 
The SKA will also provide the sensitivity to test 
whether "radio-quiet" objects like magnetars are 
in fact very weak radio emitters. The pair cascade 
processes that seem to be relevant to create the 
radio emitting plasma may be prevented if the 
magnetic field strength at the polar cap exceeds 
the critical magnetic field 

B q = m 2 e c 3 /eh = 4.4 x 10 13 G. (6) 

For magnetic fields of such strengths, other pro- 
cesses may compete with magnetic photon split- 
ting, changing the opacity of the gap region. Such 
arguments have been used to explain the absence 
of radio emission from magnetars but these 
models have been challenged by the discovery of 
high magnetic field radio pulsars |8I52| . The SKA 
can both establish whether magnetars are weak 
radio emitters as well as discover more radio-loud 
pulsars with spin-parameters similar to those of 
magnetars, so that the question about the re- 
lationship between pulsars and these prominent 
high-energy sources can be settled. As the Galac- 
tic menagerie of pulsars grows, additional classes 
of objects may also be identified, such as young 
NSs with anemic magnetic fields that appear to 
be lacking on the left side of the P — P diagram. 

The Radio- High- Energy Connection: Regularly 
pulsed high-energy emission is known for an in- 
creasing sample of radio pulsars |76| . As for the 
radio emission, the origin of the high-energy emis- 
sion is still a matter of significant debate. Likely 
emission processes for the observed non-thermal 
optical, X-ray and 7-ray emission are synchrotron 
emission, curvature emission and inverse Comp- 
ton scattering. The location of this emission 
is thought to be either in the polar cap region 
(e.g. [3S|) or in the outer gaps (e.g. ^O]). In 
both families of models, the high-energy emission 
beam appears to be much wider than the radio 
beam. With upcoming missions like GLAST, the 
high-energy emission of pulsars will be much bet- 
ter understood by the time observations with the 
SKA can be made. Then, however, the SKA can 
be used to probe the relationship between high- 
energy and radio emission, so that a coherent 
model of the active magnetosphere can be de- 



veloped. Current results suggest that the high 
energy emission may indeed be related to some of 
the radio emission properties (e.g. [62). 

A radio-high-energy link is provided by the 
phenomenon of giant radio pulses. All pulsars 
with detected giant-pulse emission are also de- 
tected at high energies |22] • Another common fea- 
ture these pulsars share is the high magnetic field 
strength at their light cylinders JT]. Interpreta- 
tion of this correlation is unclear at present. Gi- 
ant pulses often occur misplaced from the normal 
radio profile but appear to be aligned with X-ray 
and/or 7-ray emission (e.g. This suggests 

a common origin of the giant pulses and high- 
energy emission and indicates that some observed 
radio emission could be in fact a by-product of the 
high-energy radiation process. This could explain 
highly unusual high radio- frequency components, 
seen to emerge at atypical pulse phases in the 
Crab profile at a few GHz [53] . Observations also 
suggest that optical pulses of the Crab pulsars 
occurring simultaneously with radio giant pulses 
appear to be somewhat brighter than other opti- 
cal pulses |6*7] , 

With its sensitivity and monitoring capabili- 
ties, the SKA can detect giant pulses simulta- 
neously at many radio frequencies, allowing us 
to determine the spectrum and their relation- 
ship to the high energy emission. This observa- 
tional mode is identical to the requirements for 
the study of the "Dynamic Radio Sky" (see sec- 
ond contribution by Cordes, this volume), as gi- 
ant pulse emission is the prototype for transient 
coherent emission. Due to their strength, giant 
pulses can be detected from sources as far away 
as the Virgo cluster and can be used to directly 
target young NS in the Milky Way and nearby 
galaxies, as discussed in M2.4I 

3.5. Resolving Pulsar Magnetospheres 

The SKA can provide the means to actually 
resolve the pulsar magnetosphere. At a typical 
distance (1 kpc), the light cylinder radius has 
an angular scale #lc ~ O^P-D^pc ^arcsec. Ac- 
tual radio emission altitudes are much smaller 
than tlc and light-travel-time arguments sug- 
gest emission sizes of 2-ns duration nano-Giant 
pulses [22 about the size of a beach ball, though 



15 



relativistic motion toward the observer enlarges 
this by a factor 7 ~ 10 3 . Conventional inter- 
ferometry has no hope of resolving the relevant 
scales. However, interstellar scintillations (ISS) 
can probe pulsar emission regions. Just as stars 
twinkle while planets do not because the critical 
angular size for atmospheric optical scintillation 
is about 1 arcsec, the critical angular size for ISS 
~ 1CP 7 arc sec. Through appropriate measure- 
ments of the diffraction pattern caused by inter- 
stellar scattering, the phenomenon may be used 
to constrain emission region sizes and locations. 
In particular, the dynamic spectrum I(u, t) shows 
structure in time and frequency that may differ 
between different pulse components. While this 
method has been applied to a few pulsars, the 
number of objects is severely limited by sensitiv- 
ity: the pulsar's flux must be detectable in a nar- 
row frequency channel (~ 1 kHz) in a short time 
(~ 10 sec). Also, because the scintillation char- 
acteristic bandwidth and time scale are strongly 
frequency dependent, flexibility in choosing the 
observation frequency and spectrometer resolu- 
tions is needed along with large sensitivity. As 
with single-pulse studies, very wide instantaneous 
frequency coverage is important for following the 
frequency evolution of ISS and intrinsic effects. 
With the SKA, similar observations will be pos- 
sible for a large number of pulsars, therefore con- 
tributing to the final solution of the pulsar prob- 
lem. Further discussion on resolving magneto- 
spheres may be found in "The Microarcsecond 
Sky and Cosmic Turbulence," (Lazio et al., this 
volume). 

3.6. Extrasolar Planets 

The first extrasolar planetary system found is 
around the pulsar PSR B1257+12 [HSj. The sys- 
tem consists of (at least) three planets, planet A 
with approximately a lunar mass, planet B with 
a mass of 4.3 ± 0.2M® , and planet C with a mass 
of 3.9 ±0.2M e Since then, a planet has 

also been found around the pulsar PSR B1620— 26 
|68I61| . in the globular cluster M4. 

While planetary systems around main- 
sequence stars had been long anticipated and 
numerous such systems have been found since, 
pulsar planetary systems were not expected. Al- 



though it is unlikely that the number of pulsar 
planetary systems will ever approach the num- 
ber of planetary systems around main-sequence 
stars, pulsar planetary systems offer unique in- 
sights. Taken together the two pulsar planetary 
systems already indicate that planets can form 
in a wide variety of environments. Already, for 
instance, the planet around PSR B1620— 26 chal- 
lenges conventional notions about the formation 
of planets. It is thought to have been acquired by 
the pulsar during a dynamical exchange within 
the globular cluster, implying that this planet 
has existed for a substantial fraction of the age of 
the globular cluster M4. This is a low-metallicity 
globular cluster, suggesting that planets can form 
in low-metallicity environments. In contrast, 
planets orbiting main-sequence stars near the 
Sun are found almost exclusively around stars 
with solar- or super-solar metallicities |29I63| . 
which has led to the belief that only stars with 
large metal contents can host planets. Similarly, 
the presence of terrestrial mass planets around 
PSR B1257+12 suggests that terrestrial planets 
will be widespread, a prediction to be tested by 
future space missions such as Kepler and the 
Terrestrial Planet Finder 

In addition to probing the details of planetary 
formation, identification of additional planetary 
systems around pulsars could provide valuable 
clues about pulsar formation and evolution and 
probe similarities and differences between plane- 
tary systems around pulsars and main-sequence 
stars. 

The existence of planets around a pulsar is 
ascertained from the times of arrival (TOA) of 
pulses over a multi-year time span. The pul- 
sar's reflex motion in its orbit about the sys- 
tem's barycentre causes the pulse TOAs to be 
delayed or advanced relative to what one would 
predict for an isolated pulsar. Thus, finding 
additional pulsar planetary systems will require 
long, high-precision timing observations. Fig- 
ure |SJ however, demonstrates that asteroid-sized 
objects and smaller can be found with SKA tim- 
ing. 
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Figure 5. Planetary companions that can be detected 
with pulsar timing observations for different values of tim- 
ing precision. With the SKA a routine timing precision 
of <lfJ.s is anticipated for millisecond pulsars. Currently 
detected planets using optical methods are shown for com- 
parison. 



3.7. Mapping the Milky Way 

The goal of a Galactic pulsar census is to de- 
tect a large fraction of the roughly 2 x 10 4 pulsars 
that are beamed toward Earth, or about an or- 
der of magnitude larger number of pulsars from 
current and near-term pulsar surveys. This cen- 
sus will necessarily determine the DMs for all 
of the pulsars as part of the discovery process, 
along with estimates for the scattering and rota- 
tion measures (SM, RM) of a sizable fraction of 
the objects in follow-up observations. We expect 
that the SKA will enable a large increase in the 
number of pulsar parallaxes, as discussed below. 

The number of sampled lines of sight will be 
sufficiently large and dense on the sky so as to per- 
mit qualitatively new models of the Galaxy to be 
developed. The current "state of the art" model 
for the Galactic distribution of the warm ionized 
medium is the NE2001 model |14I15| . This model 
describes the Galaxy largely in terms of a small 
number of smooth, large-scale components. How- 
ever, it also incorporates a number of electron 
density enhancements ( "clumps" ) and voids along 
the line of sight to a limited number of objects 
(primarily pulsars). These clumps and voids are 



required because the various observables for these 
objects (typically pulsar DMs) cannot be repro- 
duced by only large-scale components. 

We envision that, with the Galactic pulsar cen- 
sus provided by the SKA, it will be possible to 
trace the Galaxy's structure, or at least its spiral 
structure, in a self-consistent manner, rather than 
by imposing it as has been done for the NE2001 
and previous models. Parallaxes will be deter- 
mined for large numbers of pulsars, providing cru- 
cial DM-indepcndent distance measurements by 
which to calibrate the local interstellar medium. 
Parallaxes will be measureable for some pulsars 
out to ~ 10 kpc distances. For many pulsars, 
HI absorption can be measured, yielding distance 
constraints and information on fine structure in 
the neutral ISM. At larger distances, there should 
be several pulsars per degree along the Galactic 
plane. Their lines of sight will have a reasonable 
probability of intersecting an H n region even over 
Galactic distances (~ 10 kpc). The methodol- 
ogy for constructing the model would be to insert 
clumps of increased electron density, represent- 
ing H II regions, along the line of sight to pulsars. 
The clumps would be inserted in a parsimonious 
manner so as to minimize the difference between 
observed and predicted DMs. Combined with the 
parallax measurements and constraints obtained 
from scattering measurements of pulsars and ex- 
tragalatic sources, it may be possible to identify 
patterns in the locations of clumps that would 
then reveal the spiral arms. Additional discus- 
sion of such modelling is in "The Microarcsecond 
Sky and Cosmic Turbulence," (Lazio et al., this 
volume). 

3.8. Mapping the Galactic Center 

As is the case for the Galaxy as a whole, de- 
tecting pulsars in the GC will allow the magnc- 
toionic medium there to be mapped. However, 
several novel aspects of the GC result in the po- 
tential for additional information to be obtained. 
Firstly, the geometrical weighting factor for scat- 
tering, W T , can change by large amounts for pul- 
sars at different distances along the line of sight. 
For instance, two pulsars in the GC separated by 
only 50 pc along the line of sight (< 1% of the 
total distance to the GC) may have pulse broad- 
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ening times that differ by 25%. 

Secondly, gas motions in the GC can be large, 
in excess of 100 km s _1 65 . In contrast to typical 
observations through the Galactic disk, for which 
pulsar velocities dominate, in the Galactic center, 
pulsar and gas velocities may be similar. If the 
pulsar and gas velocity vectors are anti-aligned, 
changes in pulsar DMs or scattering parameters 
could happen on much shorter time scales than 
are obtained in the Galactic disk. In 1 yr the line 
of sight through a gas cloud could sweep out of 
order 100 AU, allowing the AU-scale structure of 
gas clouds in the GC to be probed. 

The distribution of pulse periods that can be 
detected toward the GC clearly will be affected 
by scattering (c.f. Figure [3J|. Nonetheless, a 
careful accounting for selection effects may allow 
the pulse period distribution to be used to esti- 
mate the past star formation rate in the GC |5()| . 
Even modest constraints on the past star forma- 
tion rate would be valuable for assessing the ex- 
tent to which the Galaxy's nucleus has undergone 
episodes of star bursts and their magnitude. 

Finally, one aspect of the SKA's Key Science 
Project on strong-field tests of gravity is to map 
out the space-time around a black hole. At a 
crude level, the distribution of period derivatives 
of pulsars around the black hole could be used to 
map out the gravitational potential on the large 
scales around the center. More exciting is the pos- 
sibility that pulsars will be found sufficiently close 
to Sgr A* that strong- field gravitational effects, 
like strong micro-lensing or large time delays, will 
be measurable |59l84l58j . In order to exploit this 
possibility for GC pulsars, the line of sight to the 
pulsar must pass within a few Einstein radii of 
Sgr A*, or within roughly O'.'l. Consequently, 
only pulsars within the central stellar cluster, or 
those behind the GC with a favourable projec- 
tion, will be useful for such strong-field tests. 

3.9. Mapping Globular Clusters 

The dynamical and stellar evolutionary his- 
tories of globular clusters are inexorably linked 
through the initial binary fraction. The discov- 
ery of tens to hundreds of pulsars in these clus- 
ters by the SKA, combined with the precise tim- 
ing it will afford, provides an excellent tool to 



study these histories. As globular clusters are 
also thought to play a vital role in galaxy evolu- 
tion, an improved understanding of their forma- 
tion and history is of importance. Recent work on 
47 Tucanae (e.g. EHj) and NGC 6752 (e.g. [22) 
has highlighted this potential, in particular when 
the knowledge obtained at radio frequencies is 
combined with results from X-ray observations 
(e.g. J3B|). 

Using the measured period derivatives of mil- 
lisecond pulsars it is possible to measure the line- 
of-sight acceleration imparted to the pulsars, as 
this is usually dominated by the gravitational 
field of the globular cluster. When combined 
with accurate positions and, possible with the ac- 
curacy of SKA timing, accurate parallaxes, this 
will lead to the best possible determination of the 
gravitational field. Via the mass density distribu- 
tion and optical observations, the population of 
non-optically active cluster members, in particu- 
lar massive black holes, can be derived. The re- 
sultant model of the gravitational potential of the 
clusters can then be used to determine the "true" 
period derivatives of the pulsars and thus their 
spin-down energies, magnetic field strengths and 
ages and also correct any measured orbital period 
derivatives. 

The more accurate measurement of the gravita- 
tional potentials in these systems will also allow 
improved modelling of the influence of binaries 
on the collapse or otherwise, as well as determi- 
nations of how many pulsars have potentially es- 
caped the cluster and thus the initial binary frac- 
tions. It will also determine what fraction of NSs 
have escaped these systems and in turn how many 
of the NSs in the Galaxy originated in globular 
clusters. 

The stars in the globular clusters in our Galaxy 
span a large range of metallicities and a complete 
census of pulsars in a large number of globular 
clusters will therefore provide vital input in de- 
termining the stellar evolutionary history of clus- 
ters and in particular how the massive end of the 
initial-mass function depends on metallicity. Dif- 
ferences in the measured dispersion measures of 
pulsars in globular clusters combined either with 
models of the potential, or parallaxes, can be used 
to probe the gas content of the clusters f|27|'). 
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The gas, originating from the winds of evolved 
stars, is expected to be stripped from the cluster, 
each time the cluster passes through the Galactic 
plane every ~ 10 s years ([7J). The amount of 
gas determined for various clusters may therefore 
provide information about the dynamics of the 
cluster as a whole. 

The study of cluster dynamics is supported by 
the determination of the proper motions of the 
large samples of pulsars. This will allow accu- 
rate determinations of the proper motion of the 
globular clusters themselves and thus, when com- 
bined with optical proper motions, will provide 
excellent ties between radio and optical reference 
frames. 

3.10. Pulsar Demographics and Core- 
Collapse Physics 

The currently known ~ 1700 pulsars represent 
just the 'tip of the iceberg' of a much larger pop- 
ulation of perhaps 10 5 — 10 6 active pulsars in the 
Galaxy. The current sample, biased by selection 
effects and sensitivity limitations, does not allow 
answers to a number of important questions that 
include: 

• How many pulsars are in the Galaxy, what 
is the birth rate, and is there evidence for 
an increased birth rate in the past? 

• Are neutron stars formed preferentially in 
spiral arms? 

• How many pulsars are in globular clusters? 

• Do the magnetic fields of isolated neutron 
stars decay? 

• What are the minimum and maximum spin 
periods for radio pulsars? 

• What is the relationship between core col- 
lapse and neutron star birth properties? 

• How are isolated millisecond pulsars pro- 
duced? 

The properties of the population of pulsars dis- 
covered in the Galactic census f ^.lfl . combined 
with the statistics of the discovered globular clus- 
ter ( M2.3fl and extragalactic pulsars f H2.4f> . will 
unlock the answers to these questions. 

The outcome of the Galactic Census and fu- 
ture high-energy studies will establish the rela- 
tionship between the different manifestations of 



NSs, i.e. active radio pulsars, "radio-quiet" spin- 
powered X-ray pulsars, and Soft-Gamma- Ray Re- 
peaters and Anomalous X-ray Pulsars known as 
magnetars (see also M3.4|) . 

Most (if not all) NSs derive from Type II, core- 
collapse supernovae. Evidence supports the idea 
that the collapse and ensuring explosions are not 
symmetric, thus imparting a translational mo- 
mentum impulse or "kick" and, perhaps, an an- 
gular momentum kick that determines the initial 
spin state of the NS. The large translational mo- 
tions of pulsars, (V) ~ 500 km s _1 represents 
firm testimony to this picture, with corrobora- 
tion from X-ray jets that are aligned with the spin 
axes in the Crab and Vela pulsars and also with 
the proper motion vectors. While collapse kicks 
seem to under ly large pulsar velocities, the de- 
tailed mechanism for producing kicks is not yet 
identified 021- Candidate rocket effects include 
neutrino emission and mass ejection during NS 
formation, which may be altered by the magnetic 
field strength that most likely grows during the 
collapse. A slower-acting electromagnetic rocket 
(the Harrison- Tademaru effect) may occur if NS 
are born spinning fast (Pinitiai ~ 1 nis) and if the 
magnetic field is off axis. In addition it is not 
known if the velocity distribution has been sam- 
pled well enough in the high-velocity component 
to determine the maximal space velocity and thus 
the degree of asymmetry. 

With the SKA, comprehensive surveys for pul- 
sars with followup astrometry (from interferome- 
try) and polarization measurements will allow a 
much better understanding of the empirical prop- 
erties of the pulsar population with regard to 
kicks and presumably to identification of the pro- 
cesses that underly them. Information will come 
from determining parallaxes and full 3-D veloci- 
ties by using the timing and VLBI capabilities of 
the SKA. Whilst it is currently impossible to de- 
termine radial velocities for pulsars (unless in or- 
bit with an optically detectable companion), this 
will become possible with the SKA for at least 
a number of pulsars since the affordable timing 
precision will allow us to detect motion on the 
sky that is non-linear in time because it results 
from the three-dimensional spatial motion of the 
pulsar projected onto the celestial sphere. 
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4. Technical Requirements 

The key observations to achieve the outlined 
science goals can be summarized as 

• Full Galactic and extragalactic censuses of 
the pulsar population by sensitive searches 
of the whole sky, the Galactic plane, and 
specific targets such as the Galactic Center, 
Globular Clusters or external galaxies. 

• Regular multi-frequency timing observa- 
tions of the >10 4 discovered pulsars closely 
spaced in time in about bi-weekly intervals. 

• Targeted observations and dedicated stud- 
ies of individual objects and/or areas of spe- 
cial interest. 

• VLBI observations of a large number of 
sources to determine their astrometric pa- 
rameters. 

We discuss the technical requirements for these 
observing programmes in more detail in the fol- 
lowing. 

4.1. Search Observations 

Yields of pulsar searches need to consider prop- 
agation effects and binary motion in addition to 
distributions of spin periods, luminosities and dis- 
tances of target populations. For both single- 
pulse and periodicity searches, algorithms ap- 
proximate matched filtering procedures on the 
signal. 

The main requirements for pulsar searches are 

1. A configuration that allows a significant 
fraction of the total SKA's collecting area 
to be used in blind surveys that sample a 
wide field (e.g. at least 1 deg 2 at 1 GHz). 

2. Frequency coverage that allows optimal de- 
tection of pulsars in periodicity and single- 
pulse searches for most Galactic and ex- 
tralagalactic lines of sight (e.g. 0.5 to 
5 GHz) and also for the highly scattered 
Galactic Center region, which requires 9 to 
15 GHz. 

3. Wide bandwidth with adequate channel- 
ization to allow dedispersion. Dispersive 
smearing across a channel is 

8.3^DMA^ ch 
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for Ar/ C h and v in MHz and GHz, respec- 
tively. DM ranges up to 1200 pc cm" 3 
in the known sample and values as large 
as several thousands are expected for some 
lines of sight. Assuming 20% fractional 
bandwidth (an underestimate at 1-2 GHz 
frequencies) , the number of channels needed 
in post-detection dedispersion is 

Wd , .^("Si-A. (8) 
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At low frequencies or high DM, dedisper- 
sion is best done coherently, requiring ac- 
cess to the pre-detection baseband signal. 
For coherent dedispersion, the only extrin- 
sic smearing effects are from interstellar 
scattering and instrumentation. The sig- 
nal processing burden is quite large for co- 
herent dedispersion, however, and will re- 
quire state-of-the-art processors to achieve 
desired band widths. 

4. Time resolution sufficient to sample the 
shortest pulse widths. 

5. Flexibility of the correlator/beam- former to 
provide the above sampling. 

We quantify the collecting area usable in pulsar 
surveys as the fraction f c in a "core" array with 
which full-FOV sampling can be performed. In 
some designs, the core array would correspond to 
those antennas whose signals can be directly cor- 
related (as opposed to being combined in a beam- 
former first). The system-equivalent flux density 
for the core array is S sys /f c = 2kT sys / f c A c = 
0.14/" 1 Jy using the nominal SKA specification 
of 2 x 10 4 m 2 K" 1 . 

Plots of SKA performance in pulsar periodicity 
searches are shown in Figure El along with those 
for existing telescopes (Arecibo, the Green Bank 
Telescope, and the Parkes 64m). 

Single-pulse searches: Searches for giant 
pulses can be made in both blind and targeted 
searches. Because they differ from periodicity 
searches only in the Fourier analysis part of the 
latter, the same requirements on numbers of 
channels, total bandwidth and overall sensitivity 
apply. 
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Figure 6. Plot of D max vs. P for surveys at 1.4 GHz 
using the Parkes 64-m antenna, Arecibo, and the 
SKA, as labelled, and the Green Bank Telescope 
(unlabelled). The four panels correspond to different 
pseudo-luminosities Lp. In each panel a band is shown 
for each telescope, the top corresponding to full gain 
and the bottom to partial gain. For the SKA, the 
partial gain is due to only f c = 0.25 of the array being 
available for pulsar searches. For the other telescopes, 
the bottom band boundary is for 0.5 of the peak gain. 
Parameters used for the SKA, Arecibo, GBT and Parkes, 
respectively, are as follows. Bandwidth: 0.4, 0.3, 0.4 
and 0.3 GHz; Integration time: 3600, 300, 900 and 2100 
sec; Number of channels: 1024, 1024, 1024 and 96; Time 
constant: 64, 64, 64 and 250 fjs; SEFD: 0.14, 3.6, 16 
and 40 Jy. In all cases, a 10cr detection threshold for 
the harmonic sum of a Fourier spectrum is assumed. 

The decline at small periods is due to smearing of 
pulses from a combination of receiver time resolution 
and propagation effects (dispersion and scattering). 
This effect is larger for the Parkes Multibeam survey 
which used wide spectrometer channels, thus allowing 
more dispersive pulse smearing. With 1024 channels, 
dispersion is largely unimportant while pulse broadening 
from scattering limits D ma x for short periods, periods. 

The results show that the SKA can detect pulsars 
at distances 2 to 10 times larger than with Arecibo and 
can probe to reasonable distances even for luminosities 
L p < 0.01 mJy kpc 2 that are currently inaccessible. For 
L p = 10 mJy kpc 2 , the SKA can reach 100 kpc for 
longer periods. Given that some pulsars have L p > 10 3 
mJy kpc 2 at 1.4 GHz, it is clear the the SKA can reach 
distances > 1 Mpc in periodicity searches. 



4.2. Targeted Observations 

In addition to regular timing observations dis- 
cussed in i !4.3l and possible VLBI measurements 
of their astrometric parameters ( H4.4|l . many pul- 
sars will also require dedicated targeted observa- 
tions to achieve the various science goals. How- 
ever, most of the requirements of targeted ob- 
servations are identical to those for search and 
timing measurements. Indeed, nearly all aspects 
of the outlined science programmes have in com- 
mon that they require multi-frequency observa- 
tions. These are necessary to study the broad- 
band emission of pulsars, to correct for propaga- 
tion effects on the pulsed signal, and for studies 
of the interstellar medium itself. It is therefore 
desirable that the SKA has constant broad fre- 
quency coverage, and at the least the SKA should 
be able to change frequency on a regular and 
rapid (~ minutes) basis. 

For some of science goals, such as the study 
of relativistic plasma physics ( H3.4I) . simultane- 
ous multi-frequency observations with full polar- 
ization capabilities are required, often over fre- 
quency ranges much greater than a few GHz. 
This simultaneity can be achieved in different 
ways with the SKA, depending on the actual de- 
sign. Assuming that SKA frontends are made 
up of a number of different classes of receiver, 
each designed to cover a different frequency range, 
simultaneous multi-frequency coverage can be 
achieved through the ability of all of the receivers 
in each of the different classes to be operated si- 
multaneously giving full sensitivity at all wave- 
lengths. Alternatively, the SKA could form sub- 
arrays of receivers from different frequencies up 
to some limited total number of receivers. Whilst 
this second option still provides the desired simul- 
taneity, it clearly results in a loss of total sensi- 
tivity. 

4.3. Timing Measurements 

Timing observations require faster sampling of 
the frequency-time plane than search observa- 
tions because the time-tagging of pulses is de- 
sired to high accuracy. At minimum, dedispersion 
must be done as completely as possible, requiring 
either coherent dedispersion of baseband signals 
sampled before detection or post-detection dedis- 
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persion with a large number of channels across 
the measured bandwidth. Sample times needed 
are <1 /j,s for millisecond pulsars, the objects of 
greatest interest. 

Timing Precision Issues: Phenomena that 
limit time-of-arrival (TOA) precision include 

(d) 

1. Radiometer noise, which yields an rms vari- 
ation 
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where S/N is the signal to noise ratio of the 
pulse peak in a sample pulse profile and W 
is the pulse width (FWHM, in time units). 

Pulse phase jitter, which occurs in the vast 
majority of pulsars as drifting sub-pulses 
and apparently random jumps accompanied 
by large amplitude variations having a mod- 
ulation index m\ = rms/mean ~ 1. Quan- 
tifying the rms jitter of a single pulses as a 
fraction /j of the pulse width W, the rms 
TOA variation is 
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where ej < 1 measures the fraction of the 
pulses in the observation time T that are 
statistically independent. Phase jitter with- 
out intensity modulation (mi = 0) produces 
a TOA error but the converse is not true. 

Scattering-induced variations include the 
systematic delay associated with pulse 
broadening (c.f. fl2.2jl . which can vary 
with epoch, and timing variations associ- 
ated with the finite number of "scintles" 
(patches of constructive interference) in the 
frequency-time plane. Closely related are 
DM-variations associated with large-scale 
electron-density fluctuations (AU to pc; [I]) 
and variations in angle of arrival caused by 



large-scale gradients. Together, this col- 
lection of effects associated with the dis- 
persion relation for cold plasma in the 
ISM yields several contributions with dif- 
ferent frequency scaling, ranging from v~ 2 

to i/- 4 - 4 Eg. 



Pulse Polarization can induce TOA errors 
from incorrect gain calibration and from in- 
strumental polarization. For a Gaussian 
pulse of width W, the induced TOA error 
from gain variation is 
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where e is the fractional gain error on the 
circularly polarized channels (Sg/g = e) and 
7Ty is the degree of circular polarization. For 
e = 7Tv = 0.1 and W = 1 ms, the error is 
4.2 /its, for example. If there is no circular 
polarization, the error vanishes. If the re- 
ceiver polarization channels are linearly po- 
larized, a similar error ensues. Instrumental 
polarization such as cross coupling in the 
feed can induce false circular polarization 
7Ty — 2t7 1 / 2 7Tl, where 7Tl is the degree of lin- 
ear polarization, which is often quite high in 
pulsars, and n is the voltage cross coupling 
coefficient. Combined with gain errors, the 
pseudo-V yields a TOA error as above. To 
achieve better than 100 ns timing precision, 
one must have e^ 1 / 2 ^ 10~ 4 W~g , for a 
pulse width of W ms ms. If the gain precision 
is 1%, this requires cross coupling 77 < 10~ 4 , 
corresponding to a net polarization purity 
of — 40dB. This may be achieved through 
good antenna and feed performance and 
post-processing calibration. Similar num- 
bers will hold for linearly polarized feeds, 
though it is generally advisable for measure- 
ments on pulsars, which are elliptically po- 
larized, to be made with circularly polarized 
feeds. 

To obtain optimised arrival times, the above 
scaling laws imply that observations with the 
highest sensitivity are needed to minimize the ra- 
diometer noise contribution and thus favour lower 
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frequencies (except for pulsars with shallow spec- 
tra). Phase jitter, whose strength (ctoAj) varies 
greatly from pulsar to pulsar, is nearly indepen- 
dent of frequency and can be minimised only with 
longer dwell times on a pulsar. Phase jitter begins 
to dominate the contribution from radiometer 
noise when S/N ~ l for a single pulse for the case 
where /j = 1. For at least one MSP, fj « 1 (2Tj. 
so that phase jitter may be less important. In 
general, however, high-sensitivity measurements 
with the SKA will yield TOAs dominated by 
phase jitter and systematic errors from polariza- 
tion and gain-calibration effects. Contributions 
from propagation effects clearly imply that high 
frequencies are favoured. The overall optimal fre- 
quency is highly pulsar dependent and thus re- 
quires that the SKA provide a wide range of fre- 
quencies, extending to at least ~ 5 GHz. More- 
over, timing measurements at multiple frequen- 
cies allow ISM effects to be partially removed. 

While precision polarization and gain calibra- 
tions are required in order to obtain the best 
TOAs, the ability to sample the frequency-time 
plane with adequate resolution to identify indi- 
vidual scintillation features is also required to cor- 
rect for some of the perturbations to TOAs. The 
effective observing time needed per pulsar to get 
a precise TOA measurement will, however, de- 
pend on the individual pulsar, i.e. as to whether 
radiometer or pulse-phase jitter contributions to 
the TOA uncertainty are dominant. 

Multiple Beaming and Multiple Fields of 
View: Given the large number of weak sources 
to be discovered among the ~ 20, 000 pulsars 
to be timed, the capability to observe multiple 
pulsars simultaneously, either through multiple 
beams within the same FOV or by using multiple 
FOVs, is clearly desirable. Along the Galactic 
plane we expect <100 sources per 1-deg 2 FOV 
while as many beams should also be sufficient to 
sample the denser populated regions. However, 
most f-deg 2 FOVs will contain only a few pul- 
sars, so that larger FOVs or several independently 
steerable FOVs would be beneficial. 

In a worst case analysis, one can estimate from 
simulations that with only four independent 1- 
deg 2 FOVs, one needs ~ 50 days of observing 
time with full-SKA sensitivity to get one TOA 



measurement for all discovered pulsars. Assum- 
ing that the main contribution to TOA errors 
for MSPs is given by radiometer noise, a high 
S/N TOA measurement for all ~ 1000 MSPs 
would still require ^10 — 20 days with full SKA 
sensitivity. Access to only reduced SKA sensi- 
tivity and limited simultaneous multi-frequency 
coverage will increase the needed observing time 
correspondingly while observations on similar or 
even more frequent intervals are clearly necessary 
to achieve some of the key science goals. How- 
ever, in reality, optimisation procedures, such as 
the choice of lower frequencies above the Galac- 
tic plane or the selection of fewer but more pre- 
cise MSPs, will reduce the demand on observing 
time. Nevertheless, given the necessity to allow 
the SKA also to be used for other key science 
programmes, a larger number of independently 
covered FOVs of the sky appears essential. 

4.4. VLBI Observations 

The scientific motivation for obtaining VLBI 
observations of pulsars is three-fold. Firstly, from 
the requirements of the strong-field gravity key- 
science project, it is clear that uncertainties in a 
pulsar's distance and motion limit the accuracy 
to which kinematic corrections can be applied to 
observables like the orbital period derivative and 
hence to which precision tests of theories of grav- 
ity can be performed (see H3.1|l . Secondly, pul- 
sar parallaxes provide key calibration constraints 
into models of the Galactic electron density. Even 
a relatively small number of parallaxes can pro- 
vide crucial constraints on the distribution of the 
interstellar medium near the Sun |78| . Finally, 
proper motions are essential for estimating pul- 
sar velocities, which in turn provide constraints 
on the supernova collapse mechanism (see H3.10fl . 

VLBI observations require finding an optimal 
frequency that balances astrometric accuracy, 
pulsar spectra, and atmospheric and ionospheric 
phase fluctuations. Pulsar spectra generally favor 
observations at lower frequencies (<1 GHz), but 
ionospheric phase fluctuations become significant 
and the lower resolution attainable (e.g., 5 mas on 
terrestrial baselines) vitiate astrometry at these 
frequencies. Higher frequencies provide higher as- 
trometric accuracy and reduced ionospheric phase 
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fluctuations at the cost of a reduced sample size 
because of the steep spectra. 

We anticipate that the SKA's sensitivity will 
provide a broad optimization for VLBI observa- 
tions in the approximate frequency range of 2— 
8 GHz. Ionospheric phase fluctuations are re- 
duced at these frequencies while tropospheric 
phase fluctuations are not yet severe. Transconti- 
nental/intercontinental baselines provide resolu- 
tions of order 3 mas, and sub-milliarcsecond as- 
trometric accuracies are being obtained already 
[Hj. A key aspect of the SKA's sensitivity will 
be its role in finding calibrators/phase reference 
sources close to target pulsars. Chatterjee et 
al. (2004) have demonstrated that astrometric ac- 
curacy depends strongly upon the separation be- 
tween the calibrator and target source. 

If 10% of the SKA's collecting area is spread 
over transcontinental/intercontinental baselines, 
it would be some 20 times more sensitive than 
current dedicated VLBI arrays. This improved 
sensitivity not only increases the sample of pul- 
sars available for astrometric purposes but also in- 
creases the density of background sources to serve 
as phase reference sources. Finally, we emphasize 
the importance of dedicated VLBI capabilities for 
pulsar astrometry. Ad-hoc VLBI arrays can be 
constructed today with collecting areas approach- 
ing 10 5 m 2 . By their very nature, though, these 
arrays are difficult to maintain for the multi-year 
durations required for astrometric projects. 

4.5. Summary 

The outlined requirements are summarized in 
Table ^ where we quote the required tempo- 
ral resolution, sensitivity and frequency coverage, 
alongside with comments on the required config- 
uration, Field-of-View (FOV) sampling and po- 
larization capabilities. 

5. Conclusions 

Given the richness and diversity of the outlined 
science case, the research on pulsars conducted 
since their discovery appears to be only the pre- 
lude to what will be possible with the Square 
Kilometer Array. As for other areas in modern 
astronomy, the SKA will revolutionise the field 
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of pulsar astrophysics. The SKA will enable not 
only the discovery of most pulsars in the Milky 
Way but also allows present-day-survey sensitiv- 
ities to pulsars in the closest Galaxies. With the 
single-pulse search techniques, it should be pos- 
sible to detect giant pulses from pulsars as dis- 
tant as the Virgo Cluster. The impressive yield 
of 10,000 to 20,000 discovered pulsars, including 
over 1,000 millisecond pulsars, effectively samples 
every possible outcome of the evolution of mas- 
sive binary stars, thereby guaranteeing the dis- 
covery of very exciting systems. We expect at 
least 100 compact rclativistic binaries, including 
the elusive pulsar-black hole systems. For tests of 
theories of gravity of general relativity such a sys- 
tem could well surpass all present and foreseeable 
competitors. It may be possible through pulsar 
timing to precisely determine the mass, spin and 
quadrupole moment of the black hole (see contri- 
bution by Kramer et al., this volume). Such mea- 
surements would go beyond what is currently pos- 
sible with binary pulsar tests of general relativity, 
perhaps even confronting Einstein's famous "no- 
hair" theorem for black holes. The structure of 
the Milky Way and its magnetic field can be stud- 
ied in far greater detail than presently possible, 
particularly with the expected increase in paral- 
lax measurements made possible by the SKA's 
superior sensitivity. In addition, studies of the 
significant numbers of extragalactic pulsars ex- 
pected to be detectable by the SKA would allow 
measurements of the intergalactic, as opposed to 
the interstellar, medium. 
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